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We derive new bounds on hidden sector gauge bosons which could produce new energy loss mech- 
anisms in supernovae, enlarging the excluded region in mass-coupling space by a significant factor 
compared to earlier estimates. Both considerations of trapping and possible decay of these parti- 
cles need to be incorporated when determining such bounds, as does scattering on both neutrons 
and protons. For masses and couplings near the region which saturates current bounds, a signif- 
icant background of such gauge bosons may also be produced due to the cumulative effects of all 
supernovae over cosmic history. 

PACS numbers: 95.35. +d 



I. INTRODUCTION 



Despite its phenomenological success, the Standard Model (SM) is known to be incomplete. Observations of neutrino 
masses and oscillations, and several naturalness issues (stability of the electroweak scale, strong CP problem) have 
motivated theoretical extensions, many of which involve the existence of new particles and forces. Extensions of 
the Standard Model by inclusion of one or more Abelian groups are familiar in a number of contexts. Additional 
U(l) groups naturally arise from the breaking of GUT groups to the Standard Model, from attempts to resolve the 
fine-tuning of the Higgs mass, or from the possible existence of new hidden sectors in nature (for a review see [HE]). 
This latter alternative has elicited significant interest recently because it allows the possibility that new light gauge 
bosons associated with the new U(l) groups may connect a hidden sector to the observed sector, allowing for the 
existence of associated dark matter particles which may have enhanced couplings because of mixing between the new 
gauge bosons and photons. 

If the boson were massless then it would mediate long-range forces, which would then necessitate a very small 
coupling to SM matter or that the new group be sequestered from the SM. Massive gauge bosons are not as strongly 
constrained, but even for additional U(l) sectors that are not charged under the SM there will still arise a coupling to 
the SM through kinetic mixing. That is, for a U(l)' field strength F'^" with an associated gauge boson A', and the 
hypercharge field strength the Lagrangian density will contain the gauge invariant renormalizable terms [3] 

C D \f^F^ + e fF Y ^F^ (1) 

If the U(1)y is embedded in a GUT theory, Planck-suppressed operators or loop-suppressed mixing from heavy split 
multiplets can induce ey ~ 1CP 8 — 1CU 2 . In the context of string theory, the possible range of ey is much larger, 
with estimated values of ey ~ 1CU 17 — 10 _s from compactifications of the heterotic string, or ey ~ 1CU 12 — 1CU 3 in 
type II scenarios [1H7|- At low energies, the kinetic mixing can be removed by considering a shift in the photon field 
—> A^ + eA ,fl , inducing e-suppressed electromagnetic interactions of the A' with strength ee: 

C D eeA'^ M , (2) 

where e = ey cos 9w, with weak mixing angle 6w 

As alluded to above, hidden sector models have garnered much attention recently because if there are dark matter 
particles in these sectors that are charged under a new dark gauge symmetry, the possibility exists for increased self 
interactions which could result, via mixing with the standard sector, in the possibility of an enhanced annihilation 
signature which was thought might explain anomalous cosmic ray data. Dark matter particles in the U(l)' sector 
could have a mass at the weak scale, with the A' mass suppressed by ^[i down to the MeV-GeV scale (there exists a 
rather larger literature on such enhancements, see for example [THIS] for models and constraints). 

Of course, if the new gauge bosons are sufficiently light, they can also produce observable signatures, using direct and 
indirect terrestrial probes. Beam dumps, for example, which are sensitive to the possibility of new light penetrating 
particles such as axions provide stringent constraints |27l I28j . In addition, one must consider the impact of new light 
gauge bosons on sensitive atomic probes such as the anomalous magnetic moment of the muon [29] . 
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Supplementing these constraints on a hidden sector are constraints from astrophysics. If hidden sector gauge 
particles mix with photons and then escape from a star, then this new energy loss mechanism can dramatically affect 
not only stellar structure but also stellar evolution. Perhaps nowhere is this more dramatic than in the case of core 
collapse supernovae. As a result, supernova cooling constraints (namely the energy loss observed from SN1987a) have 
long been employed in order to constrain particle couplings and masses |30H33j . 

Early work was concerned principally with constraints on nearly massless axions emitted in nucleon-nucleon 
bremsstrahlung processes inside of the hot (T ~ 30MeV) neutron star born from an associated supernova. One 
can adapt this line of inquiry to constrain the masses and mixings of a hidden sector U(l)' gauge boson, A', which 
may couple to Standard Model-sector charged particles via the previously mentioned kinetic mixing. An estimate of 
the energy loss of such a process was given in [27], but, as is mentioned in that work, this estimate was expected to 
be accurate at best at the order of magnitude level. 

It is the purpose of this work to provide a more thorough accounting of this process and thus derive more accurate 
constraints on the masses of such particles, m,A' , and their couplings, e. Because models of interest for dark matter 
involve gauge boson masses which exceed 1 MeV, we restrict ourselves to this range here. Moreover, for much lighter 
masses, other considerations arise. In particular, matter effects as particles traverse the high density collapsing core 
can be large enough to produce oscillation phenomena reminiscent of that for neutrinos which can cause enhanced 
trapping of such particles and thus may obviate the constraints we derive on heavier gauge bosons |34j . Dark matter 
particles in the U{1)' sector are generically heavy and would not contribute to the cooling of the supernova. In 
scenarios where dark matter particles have masses below ~ 10 MeV, that could explain the 511 keV line emission 
from the galactic bulge HS1 , additional cooling channels could be present [Ml GH] > which we do not consider 
below. 



II. FREE STREAMING 



We are interested in exclusion regions in the e — vtla 1 plane. What one expects is that there will be an upper bound 
on e above which the vector particle production will exceed the bounds from SN1987a, extending up to a value when 
the vector particles would be strongly enough coupled to be trapped inside the supernova core. This is what we will 
call the free streaming region. In the first sub-section we will outline calculations for the process p + p — > p + p + A' , 
and then extend this to include the p + n — > p + n + A' process as well. 



A. The proton only case 



The central quantity to calculate is the emission rate of the vector particles, which is labeled Qij where the 
subscripts stand for the possible nucleons participating in the emission process. In the hot medium with temperature 
of a few MeV and densities typical of the accretion disk, the nucleons are nondegenerate and nonrelativistic. The 
main emission process is the nucleon-nucleon-dark boson bremsstrahlung. In studies of axion emission one usually 
neglects the mass of the emitted particle. Since our dark gauge boson can have a mass comparable to the temperature 
of the SN medium, this approximation is not justified and we follow the kinematical analysis of |38j . The details of 
this somewhat lengthy calculation in the one-pion-exchange (OPE) approximation are given in the appendices. If we 



take the expression for Q pp given in Eq.(B.23l, and use the fact that the A' couples to protons and charged pions 
with strength = g p = ee, we can write: 

n 2 T 25 n 2 

Qvp= ^ 5 ^ 5 2f pp e 2 e 2 X k (y, q ), (3) 

where Ik denotes the phase space integration in terms of the dimensionless variables y — m^/mpfT, and q = rriA' /T. 
We have also replaced the baryon number density, Ub, by its mass density through = pjm^. For the pion- 
nucleon coupling we take f pp ~ 1, resulting in w 15 in Eq. pi). In our calculations we will use the typical values 
p = 3 x 10 14 g/cm 3 = 3 x P14 and T ~ 30 MeV for the density and temperature in the interior of a supernova. We 
then constrain the parameter space by imposing that the luminosity due to vector emission must be less than: 

C v < 10 53 erg/s « 4.1 x 10 37 MeV 2 , (4) 

which is roughly the energy lost in neutrinos |39j . 

Since our expression for Q pp , Eq. ([3]), is the luminosity per unit volume, we need to integrate it over the volume 
where the vectors are emitted. We then assume that T ~ 30 MeV holds within a central spherical region of 1 km in 
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radius (if we take the SN to have a mass ~ 1.5M Q , then it would have a total radius of ~ 13 km assuming the same 
constant density). 

Putting all this together, we find that the luminosity due to vector emission 

Q PP V = 1.05 x 10 48 p\A\ v r\ m e 2 X k {y, q) MeV 2 « 4.7 x 10 52 e 2 Z fe (y, q) MeV 2 , (5) 

cannot be larger than C v in Eq. Q. 

So our constraint on the coupling can be written as: 



3.9 x 10- 11 /8.8 x 10~ 16 



and we then generate an exclusion bound with y fixed, while varying q for different masses ttia 1 ■ This will give an 
upper bound on the coupling e in terms of tcia 1 ■ 

B. Including neutron processes 

Next we will include processes with both protons and neutrons. The expression for emission due to the process 



p + n^-p + n + A' is given in Eq. (D.l ). This situation is slightly more complicated as we now have five integrals 



contributing in the phase space integration. We can use Eq. (16| to find the limits, if we add a multiplicative factor of 
8 in Q. This comes about as follows: there is a factor of 4 = (\/2) 4 from isospin requirements that the neutron-proton 
coupling to a charged pion is f pn — y/2f pp ; there is no longer a symmetry factor of 1/4, since proton and neutron are 
not indistinguishable; finally, we had a factor of 2 from y 2 + gi in front of 1%, but now yg = 0. The results are shown 
inFig.[T} 

III. INCLUDING DECAY 

Naively, larger couplings that do not satisfy Eq. ^ would be excluded. But, the free streaming limit is not the only 
region in parameter space where constraints will arise. In general, due to the couplings in Eq. pi), the dark gauge 
boson decays back into leptons and other SM particles within a distance: 



3E A > 



N eff m A ,ae 2 



(7) 



In order to implement this, we will follow the strategy outlined in [37] and modulate the emission amplitude with an 
exponential damping factor, 

e -10fcm/i _ (g) 

The assumption is that the products of the decay remain within the SN core and do not contribute to the cooling [33] . 
As the mixing parameter e increases, the decay happens earlier and the gauge bosons are not effective at cooling, as 
seen m Fig. [1] Therefo re we would expect that above some value of e, the emission constraints will no longer hold, 
and we will find a lower bound on the coupling, creating a region bounded above and below once decay is included 
with free streaming. 

In terms of rfe m and the decay length, Eq. ([T]), the exponential parameter can be written as: 

C - 2 .4X 10 "^^ (g) 

where x = Ea> /T. 

This exponential factor enters the phase space integration, and as a result, the function I(y, q) becomes a function 
I rfec (2/ ) q J e ) rfc m ), but we can still use Eq. M with l dec . When generating the plots, we took N e ff = 1, which 
corresponds to the situation where only the aecays A' ~ > e + e~ are possible. We also take = 10 for the typical 
distance where decay products are trapped, and, as in the previous section, we assume that vector particles are 
produced within the inner km. 
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FIG. 1: The total excluded region, adding contributions from both the pp and pn processes 



IV. TRAPPING LIMIT 



Thus far we have examined the energy loss due to vectors which are emitted and subsequently decay outside of the 
supernova. However, as is well known, above a certain coupling, these particles can be trapped within the supernova. 
This will create a lower bound on the coupling e, above which the trapping will not allow vector emission to contribute 
to cooling the supernova. 

If the new particles generated in the nucleon interaction processes subsequently interact strongly enough they will 
thermalize, and will be emitted from a spherical shell where the optical depth is roughly r rs 2/3. If the temperature 
is T x , the luminosity will be given by the Steffan-Boltzmann law: 



Here a is the Steffan-Boltzmann constant, which is 7r 2 /60 for photons and 

2 

° = 120' 



(10) 



(11) 



for a new particle with g degrees of freedom. 

In principle, one should determine r and T x , but following [315] we assume r ~ 10 km. This is a good approximation, 
because the density of the protoneutron star falls abruptly near that radius. 

The requirement that the luminosity in new particle emission be less than 



L < 10 w crg/s « 4.1 x 10 d 'MeV 



(12) 



translates into a bounded relation between the temperature and the coupling 

T x < 11.4 g- 1/4 MeV. (13) 
To calculate T x one needs a model for the temperature and density above the settled inner SN core, for which we 
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turn to [3T] : 



P(r) = PR 
T(r) = T R 



R 

r 

p{r) 
PR 



1/3 



(14) 



where pu = 10 14 g/cm 3 , R — 10 km, Tr = 10 — 20 MeV, and n ~ 5 is a relatively large number. 

As shown below, given a model, the opacity n(p, T) can be expressed in terms of r. It is then straightforward to 
compute the optical depth, 



Kpdr. 



(15) 



One then solves for r x from r (r x ) ss 2/3, and uses Eq. (14 1 to find T x . The temperature just found is used in 
Eq. (10) or compared to Eq. (13). 

To find the opacity as a function of p and T, one starts from the reduced Rosseland mean opacity: 

1 1 



K x p 4aT 3 
15g 



3tt 4 T 5 



dE A ,l EA ,/3 EA ,d T B E 
E\, (E\, 



dE A ,l 



E.r 



m A>) z Ea ' /T 



(e E A>/T _ i) 



where we have used a = 7r 2 /15, the velocity Pe a , = \/l — ijnA' I 'Ea 1 ) 2 \ and the definition of B Ea , is given by 

g E%(E%-m 2 )V 2 



Br 



2vr 2 



oE A ,/T _ i 



(16) 



(17) 



Since the production of bosons involves a Bose stimulation factor, we need to add a factor of (l — q~ e a>/ t ^ under 
the integral above. This gives the reduced opacity (often denoted as ft*), which is the one that we have been using. 
Dropping the superscript *, we find that: 



1 150 

K x p 8ir 4 



dxl 



x 2 (x 2 — q 2 ) e 2x 



(18) 



which reduces to the expression for axions |39j for g = 1, and q = m x /T = 0. 

In the corresponding calculation for a massless axion, the mean free path can be obtained as follows: starting from 
the energy-loss rate, Q x , one removes the phase-space integral J °° dE a> T7t E\, / '(27r) 3 , and a factor of Ea', since Q x is 
an energy loss rate. Additionally a factor of e Ejl ' l T must be included to account for the detailed-balance relationship, 
and this gives , ■ 

In our case, we have to take into account that the boson under consideration is massive. Therefore, the vector 
boson phase-space integral, is J d 3 p — j <inp 2 dp, with p — J E 2 A , — m 2 A , , and Airp 2 dp — 4ny/E A , — m 2 A , E A' dE a* ■ 
Hence, the inverse mean free path is obtained by removing from the luminosity, Qij, the factor: 



dE A , 



AirE A < \[E 



(2tt) s 



-^E a - = T 4 



dx 



4irx 2 y / x 2 — q' 2 
(2^)3 



(19) 



and adding the detailed balance term. 

Our expression for Qij is once again given by Eq. (B.23) for the p + p^p + p + A' case, and in Eq. (D.l ) for the 
p + n^-p + n + A' bremsstrahlung. We can write them both as follows: 



Qij — 



2 T 2.5 2 



32 7 r 1 - 5 m^ 5 



dudvdx\/uve 



-8(u — V — x)l; L . 



(20) 



where are functions defined in the appendix. The factor rj^ is equal to 2 for pp scattering, and rj pn = 16, with the 
usual isospin and non-identical particles enhancement in pn scattering. 
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Next, we perform the u- integration with the help of the delta function, and we remove a factor of Ea> = Tx from 
the x-phase-space integration. Additionally we must include the detailed balance term, which cancels an e~ x that 
appeared after the delta function was used. 

Once all of these machinations are completed, the mean free path corresponding to each channel is found to be 

! _ ^-^.5^2 (2?r) 3 



x 3 



1 f°° 

- / dvy/(v + x)ve- v I l3 . (21) 



=Sij{x;q,y) 



We can now plug this expression into Eq. (18), and manipulate it to find 

P 



16m 5 N 5 13 Ti-B : 



(22) 



where we defined the dimensionless opacity: 



K V ~ JW4*,.. / dX 1A3 c t-\> K 26 ) 



S^VijJq (e x -lf S^xy 

and S is defined above. We find the total contribution from both channels using the fact that inverse opacities add: 



One then defines the quantity tr = krprR. Substituting pr and Tr in Eq. (22), we can find tr, and 



\f ■ (24) 



Finally the opacity is found from Eq. ( 15 ) to be 



/ rji \ 9/2-3/n 

tr J-x 



2' 



This expression can now be used to bound the coupling e as a function of the dark gauge boson mass rriA> by 
requiring t x < 2/3. For the case considered in this section, of trapping without including decay, we find numerically 
the exclusion region shown in Fig. [T] 

In a more realistic scenario one must take into account both trapping as well as decay. This will lead to an exclusion 
region which is the intersection of the regions from trapping and decay alone, as either process will ensure that no 
additional energy loss will occur. For small masses and stronger couplings, trapping constrains the excluded region of 
parameter space. Conversely, for larger masses and smaller couplings, decay becomes most important. The combined 
region is shown in Fig. [l] 

In Fig. [2] we display 07] the combined region along with constraints due to other sources such as beam dump 
experiments E774, E141, E137, contributions to the anomalous magnetic moment of the muon and the electron, 
and a e , and BABAR bounds from upsilon T(3S I ) decays. 



V. A DIFFUSE DARK GAUGE BOSON BACKGROUND? 



With an average event rate of approximately 10~ 2 yr" 1 per galaxy, there have been over 10 18 supernova explosions 
over cosmic history within our observable horizon. As was recognized even before SN1987a allowed more accurate 
estimate of neutrino production rates in supernova, and before recent supernova surveys allowed a more careful 
determination of event rates, the cumulative effect of such supernovae can produce a significant diffuse background of 
particles @D]. In particular a neutrino background in the MeV range with a flux in the range of ~ 10 2 cm _2 sec _1 has 
been predicted. If a new dark gauge boson exists which is not ruled out by the arguments presented here, but which 
nevertheless contributes significantly to supernova cooling, a similar background of such particles would be expected. 
Indeed, similar considerations have recently been recognized for a possible diffuse axion background from supernovae 
[41] . Because these particles are massive, the dominant contribution to energy loss will come from their mass as long 
as it exceeds the mean temperature near the surface of the supernova neutrinosphere. As a result we would expect a 
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FIG. 2: The total excluded region for MeV to GeV dark gauge boson masses is shown including the supernova bounds calculated 
in this work, shown in blue, along with constraints due to other sources such as beam dump experiments E774, E141, E137, 
contributions to the anomalous magnetic moment of the muon and the electron, a M and a e , and BABAR bounds from upsilon 
T(3S) decays. 



diffuse background today of dark gauge particles of mass rriA> as large as w 10 2 (MeV/mA') cm _2 sec~ 1 . Their net 
contribution to the overall mass density of the universe would be negligible, but depending upon their interactions, 
such a background might be detectable. Moreover, if these particles are not stable, but decay into Standard Model 
particles, they would contribute to a diffuse cosmic ray background. We are therefore currently examining whether 
various direct and indirect constraints may be relevant for the possible detection of such a background [T^] . 



VI. DISCUSSION 



The calculations described here produce more accurate, and more importantly, significantly enhanced constraints 
on light gauge bosons that may arise from hidden sector models which predict novel types of dark matter. Mixing 
with photons produces results in bremsstrahlung production by scattering of both protons and neutrons, and a careful 
consideration of the resulting constraints from observations of SN1987a, involving both considerations of trapping and 
possible decay allow us to rule out a significantly larger region of masses and mixings than was previously estimated 
)27j . extending a constraint on an effective coupling of less than 10 -10 for a significant range of gauge boson masses. 
We also note that for masses and mixings near those that saturate the supernova cooling bounds, a significant diffuse 
background of dark gauge bosons will be produced by the cumulative effects of all supernovae over cosmic history. 
Whether such a background might be detectable is of some interest, and is the subject of ongoing investigations. 
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Appendix A: Amplitude for the process p + p^-p + p + A' 

Nucleon- nucleoli bremsstrahlung is the primary mechanism for dark gauge boson emission from the inner core of 
the SN. Like in the corresponding axion calculation, we calculate the rate for this process in the one-pion-exchange 
(OPE) approximation, which is sufficiently reliable for our present purpose @3] (although it could overestimate the 
emission rates by a small factor at lower temperatures |44|). Since the hidden sector boson only couples at tree level 
to SM particles with electric charge, we need only consider nucleon scattering processes which involve protons. We 
have eight diagrams for the OPE proton scattering processes with vector bremsstrahlung, p + p — > p + p + A' and 
five diagrams for the p + n — > p + n + A' process, where A' is the hidden sector U(l)' gauge boson with mass mA>, 
and we will take all nucleon masses to be m^r. In this Appendix and in Appendix B we will examine the process 
p + p^-p + p + A', with a treatment of the process p + n — > p + n + A 1 given beginning in Appendix C. 

Following [32] these diagrams, as shown in Fig. ([3]), are designated a, b, c, d, a',b',d,d', which gives a total of sixty- 
four terms from the squared matrix element 

^MM f = Y J \Ma + M b + M c + M d -M' a -M' b -M' c -M' d \ 2 (A.l) 

5 S 

The plus sign is for the direct t— channel diagrams and the minus sign is for the exchange u-channel diagrams. As 
required by gauge invariance, the amplitude satisfies 

k%M^ = 0, (A.2) 

where M^ represents the Feynmann amplitude without the external polarization vector. 
The matrix elements for each diagram are given by 

u{p i )^u{p 2 )u{pz)i{f ?) + $A> + m N)l5u(pi) 
u{p3)lbu{p2)u{pA)i{f i + $A' + ^Af)75"(Pl) 
w(P3)75-"(piMp 4 )^ 4 + %A' + m N)lsu{p 2 ) 

u(p4h5u(pi)u(p 3 )^(p' 3 + $a' + m N )j 5 u(p 2 ) 
-u{p4)l^u{p 2 )u{p z )-fr {f i - $ A , + m N ){u{px) 
u{p3)lbu{p 2 )u{p i )-i b (f 1 - $A' + mN){u(pi) 
u{pz)lhu{pi)u{p i )^ b {f 2 - % A , + m N ){u(p 2 ) 
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where f pp is the coupling for the p—p—ir vertex, and we have used the momenta definitions fc = p 2 —p4 and I = p 2 —p3- 
With these definitions we find the kinematic relations 

2 I k m A f 
P\-P2 = -K-1 + m N - 
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FIG. 3: Diagrams for the process p + p — > p + p + A' . Following [32] , these are labelled from top left to bottom left as a, b, c, d, 
and from top right to bottom right as a' ,b' ,c' ,d! . 

As in previous work with massive axions |38j , we will work in the approximation that the nucleon mass is much larger 
than both the temperature and the vector mass, and that the direct and exchange four-momenta transferred is larger 
than the vector mass and momenta, k 2 , I 2 ^> k\, , m 2 AI , E\, . This is due to the fact that the kinetic energy of the 
vector particles should roughly be a factor of a few times the temperature. We can then use the simplifying relation 
Pi ■ kA> — vtlj^Ea 1 ■ Finally, we use the approximation m 2 N 3> |k| 2 , since |k| 2 ~ 3m nT <^ m 2 N . 



Squaring the matrix elements, we find for the process p + p — > p + p + A' 



\M\ 2 P+P = 



' N 



64k 2 



C fc |k|' 



C H (|k| 2 |l| 2 -2|k-l| 2 ) 
m% E\, V(k 2 + m 2 )(k 2 + m 2 ) ' (l 2 + m 2 )(1 2 + to 2 ) (k 2 + to 2 )(1 2 + m 2 ) 



(A.3) 



where C k = f^ p {g 2 a + gf), Q = fp P (gl + gj - 2g a gf}), and C u = f^ p {g 2 a +g}- 2g a gp)- Here the g, are the couplings 
between the A' and the i'th particle; which in this case is g a = gp = g p = ee. Thus, Ck will be the only non-zero 
contribution. 



Appendix B: Phase space integrations for the process p + p — > p + p + A' 

The energy emission rate is given by 

Qa = JdKSY; \M\ 2 (2ir) 4 E A ,6 4 ( Pl + p 2 - p 3 - p A - k A ,)fJ 2 , 

spins 

where 



dn 



d 3 p 1 d 3 p 2 d 3 p 3 d 3 p4 d 3 k A > 
(2n) 3 2E 1 (2tt) 3 2E 2 {2ir) 3 2E 3 (2n) 3 2E 4 (2tt) 3 2E a , 



(B.l) 



(B.2) 
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The initial-state nucleon occupation numbers are given by the non-relativistic Maxwell-Boltzmann distribution 



tib ( 2tt \ 



3/2 



-p 2 /2m N T 



(B.3) 



and S is the symmetry factor, 5 = 1/4 for the pp process and S = 1 for the np process. In the non-relativistic 
approximation, Ei ~ m^. 

Following [32, 38, 45 we go to the ccntcr-of-momcntum coordinates 



Pi 

P2 
P3 
P4 



= P+Pi 

= P-Pi 

= P'+P/ 

= P' - P/- 



(B.4) 
(B.5) 
(B.6) 
(B.7) 



In the non-relativistic limit, a typical nucleon with kinetic energy E^ in has momentum |p| ~ ^/2m^E^ in -Et^ 
\&.A' ■ Then, the three dimensional delta function will enforce P = P'. 
One defines the dimensionless variables 



m N T 



Efr_ 
T 



y 



7T „ 



m,A> 
T 



(B. 



as well as z = cos(#jf), where Oif is the angle between and pf. The velocity of the emitted boson, £, can be written 

as 



c 2 = 



k^ 2 

Ej, 



= (1"^)- 



(B.9) 



The momenta in Eq. (A.3) can then be written as 



Ikl 2 


= toatT(u + v 


— 2z\/uv) 




(B.10) 


111 2 


= toatT(u + v 


+ 2z^/uv) 




CB.H) 


k| 2 + ml 


= m^T(u + v 


— 2zyuv -+ 


-y) 


(B.12) 


111 2 +"4 


= m^T{u + v 


+ 2z\Juv 4 


-v) 


(B.13) 


|k-l| 2 


= (m N T) 2 (u - 


-vf. 




(B.14) 



The delta function in the emission rate, Eq. (B.ll, becomes 

$ (pi +P2 - Pa -Pi - k A >) = S (3) (Pi + P2 - P3 - Pi - k A >)5(u-v - x)/T, 



(B.15) 



and we use the three-dimensional piece to do the P4 integration. 

We now change coordinate bases from pi, p 2 , P3 to the Pi, p/, P system, including a factor of 8 (the inverse 
Jacobian). 



Using 



we obtain 



d p j; e 



J I J 2 



-p 2 /m N T _ 



2tt y 



4 \tojvT / 



e -« e -P /™«T 



47T 



p 2 e -p?/ m „T rfpi = 27T (m N T) 3 ^ 



/ue u du 



J d 3 p f = 2ir J p}dp f J dz = TT(m N T) 3/2 J ^/vdv J dz 

d 3p c -P 2 /m«T = ^ f p2 e -P 2 /m N T dp = ^ mNT ^/2V^ 



(B.16) 
(B.17) 

(B.18) 
(B.19) 



Gathering the factors of E A i in the initial emission expression, Eq. (B.ll, the denominator from the phase space 
integration, and the amplitude squared, Eq. (A.3), we are left with 



/ 



d 3 k A , 



1 



(B.20) 
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Using 



k A , 2 = C?E\, 



(B.21) 



we have 



d 3 k A , 



E 2 
A' 



4ir / dk A > 



E\, 



4tt / dk A ,( 2 = 4ttT / dx 



\J x 2 — q 2 



(B.22) 



Using the above relations for the phase space integration factors, we insert \M.\ 2 given in Eq.(A.3l into Eq.(B.l|, 
and find that the emission rate for the p + p — > p + p + A' process is given by 



a 2 j>2.5 n 2 

P^- / dudvdzdxWuve 



ipp 32^^ 



8{u-v- x) x (I k + Ii+ I k i + I k4 ) 



(B.23) 



where a T = (2mjv/ f PP ) 2 /47T and the various integrand pieces are defined as 



(u + V — 2z^/uv + y) 2 

(u + V + 2Zy/uv) 2 



:(u + v — 2zy/uv) 



h = f 4 (fjl+90-^9 a 9p)- f ,, — ,,■ 

' (u + v + 2zy/uv + yy 

»a, 9 9 n ( u + f) 2 — Az 2 uv) , „ , — . 

1 (u + v + y)^ — 4:Z z uv 



-2/ 4 (<? Q + flg - 2g ag/3 ) , \_ ; 2 (tt + « - 2z^v) 
' (u + v + y) 2 — Az z uv 



(B.24) 



Appendix C: Amplitude for the process p + n^p + n + A' 

The requisite matrix elements for the process p + n — > p -\- n-\- A' are 

Ma = |kp+^2m^(^) *t^(^^+%)7 5 »i^ (CI) 
M c = - |k|2 '^ 2 J^ t (^f) fc&^fc " ^ +m.)/W- 2 (C2) 

1 1 / 2772 \ ^ 

^ = lip + m^m^i^ J fr&*»Wa + t A > + m N )^u^m (C.3) 
M c , = m% 2m ^ A/ {^)g P f 2 pn u^{ tl - U + m^fisT 5 ^ (C4) 

1 1 / 2771 \ ^ 

Me ' = + (Z-a)2- m j(-^J 9,f 2 pn ua b u 1 u^u 2 {k A ,-2l)-e (C.5) 

where the coupling f pn is the coupling for the p — n — it vertex, which is related to f pp as f pn = V2f pp by isospin 
invariance. As in the pp case, gi is the coupling between the A' and the ?'th particle. For this process we have 
9n = g P = ee. 

The products of matrix elements for diagrams with bremsstrahlung originating off of external legs will produce 
terms in \A4\ 2 identical to those found in Appendix A for the p + p^p + p + A' process. However, we see that we will 
also find matrix element products which include the factor Me' which arises from internal bremsstrahlung. These 
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new combinations are given by 

Me>M\, = -a e ' e 'Tr[(/ 4 + m A r) 7 5 (|) 1 +m A r)7 5 ]Tr[(^ 2 + m A r) 7 5 (^ 3 + m 7 v)7 5 ] 

x (k A > - 21) ■ (k A > - 20 (C.6) 

a ,, = 1 1 2 f 4 fC7) 

(|l| 2 +m 2 ) 2 {{ i _ kA ,y _ m 2y m 4 ^Jm K ■ 

M a M\, - -Ctae, (C.8) 

x Tr[(f 3 + m N ){^ A , - 2/)(jf 3 + + m N ) 1 b {f 1 + m N ) 1 b (f i + m N )^(f 2 + m N )j 5 } 

1 1 1 8m% 4 

aae ' {W+ml){\^+ml){{l-k A ,y-ml)E A ,mi 9 * 9pJ v n 

M c Ml = - a<x , (CIO) 
x Tr[{f 3 + m N ) 1 5 (p i - % A , + m N )($ A , - 2/)^ + mjv)7 5 Q* 4 + ™jv)7 5 (?> 2 + ™w)7 5 ] 

(C.ll) 

1 1 1 8m% 4 

(|iP+m2)(|k|2 + m2)((/-fc A ,) 2 -™2)i? A ,m4 5 ^" 

Mb'Me' = -av e ,Tr[(f 3 + m N )(# A , -2/)(j> 3 + j^, + m A r) 7 5 (^ 2 + m A r) 7 5 ] (C.13) 

x Tr[(|i 4 + toat) 7 5 (|) 1 + toat)7 5 )] 

1 1 8m 3 N 4 rri^ 

ab ' e ' - (|ip + m2)2(a-^) 2 -^)^m4 ff ^« ^- i4j 

Mtf-M^ - -ace' Tr[(^ 4 + TOAr) 7 5 (|» 1 -jlS jl /+m JV )(^ -2/) (/ 1 +to J v) 7 5 ] (C.15) 
x Tr[(f 3 + m N )j 5 (f 2 + m N )j 5 )} 

1 1 8m% f i (r 1fi x 

ac ' e ' " (|l| 2 +m2)2(a-^) 2 -^)^m4 5 ^« ^- ibj 

(C.17) 



Putting everything together we find 

64m 2 r |k| 2 / C fc |k| 4 , 



2^ 'np-ynpa ^2^4 ^(| k |2 +m 2)2 + 



|2 +m 2)2 



128m^|k-l|/ C^(|k| 2 |l| 2 -2|k-l| 2 ) 



+ 



EA'fn% \ (|k| 2 + m 2 )(|1| 2 + m 2 )((Z - M 2 - ™ 2 ) 



2CWIH 4 



2 +™ 2 ) 2 ((Z-fc A 2 -^) 
64m^ C E >E'tt\ 6 



m% (|l| 2 + m 2 ) 2 (a-fc^) 2 -^) 2 

where 



(C.18) 



Cfc 


= f 4 o 2 = f 4 e 2 e 2 


(C.19) 


c, 


- « = 


(C.20) 




<-4 2 2 
— 9n9pJpn ~ Jpn e e 


(C.21) 


Cb'e' 


i-4 j-4 2 2 
— 9n9pJpn ~ Ipn e € 


(C.22) 


Ce'e' 


= Q 2 f 4 = f 4 e 2 e 2 

-yirJpn J pn^ 


(C.23) 
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FIG. 4: Diagrams for the process p + n — >• p + n + A' . From top left to bottom left, the external bremsstrahlung diagrams are 
labelled a and c, while from top right to bottom right, the external leg bremsstrahlung diagrams are labelled 6' and c' . The 
internal bremsstrahlung process is labelled e'. 



Appendix D: Phase Space Integrations for n+p — ¥ n + p + A' 

We can now compute the phase space integrations in the same manner as in the p + p — > p + p + A' case. We find 



Q 



pn 



J dudvdzdxyfwvc u y 1 — -^8(u — v — x) 



X {Ik +Ii + Iae'h + lAE'k-l + Ib'E' + Ie'e 1 ) 
where (note the x dependence of the last four) 

(u + v — 2zy/uv) 3 



h — Ck 



II 

Iae'h 

lAE'k-l 



= Cl 



(u + v — 2z^/uv + y) 2 
(u + v + 2zy / uv) 2 



,„ (u + V — 2zy/uv) 

(u + v + 2zy/uv + y) 2 

((u + v) 2 — 4z 2 uv) (u — v) 

((u + v + y) 2 — Az 2 uv) (u + v + y + 2zy/uv) 
(u-v) 3 1 



2xC A e> 
-AxCae 



((u + v + y) 2 — 4z 2 uv) (u + v + y + 2z^pwv 



(u + v + 2zy/^) 2 
Ib'E' = 4xC b >e> 7—, - , n _ 7= , . ^J u-v) 



Ie'E' = x 2 Ce'E' 



(u + v + 2zy/mi + y) 3 
(u + v + 2zy/uv) 3 



(u + v + 2zy/uv + y) 4 



(D.l) 
(D.2) 

(D.3) 
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